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FOREWORD 

This report describes a portion of the results obtained on NASA Grant NSG 
3044. This work was done under subcontract to the University of Illinois » 
Urbana, with Frof. S.S* Wang as the Principal Investigator. The printe gran*- 
tee was the Massachusetts Institute of Technology! with Prof. F.J. McGarry as 
the Principal Investigator an^ Dr. J.F. Mandell as a major participant. The 
NASA - LeRC Project Manager was Dr. C.C. Chamis. 

Efforts in this project are primarily directed towards the development for 
finite element analyses for the study of flaw growth and fracture of fiber 
composites. This report presents a convenient method of analysis for cracks 
emanating from holes In anisotropic materials. 
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ABSTRACT 

A mechod of analysis has been developed for craeks emanating from a cir- 
cular hole In an off-axis unidirectional f Ibcr-relnforced composite. This 
analysis method Is formulated on the basis of conservation laws of elasticity 
and of fundamental relationships In anisotropic fracture mechanics. The prob- 
lem is eventually reduced to a system of linear algebraic equations in mixed- 
mode stress Intensity factors. Superiority of this current method to other 
approaches in investigating the problem with a very complicated crack geometry 
and material anisotropy Is demonstrated when used in conjunction with any 
numerical method such as a finite element analysis. This method Is found 
to be convenient and accurate. Mixed-mode stress intensity factors and the 
associated energy release rates In the crack problem arc determined for the 
composites with various fiber orientations. Solutions for both single and 
double cracks emanating from the edge of a hole in the composites are pre- 
sented also to illustrate the fundamental nature of the problem. 



1. INTRODUCTION 


Fracture of fiber *-reln forced composite materials containing stress con*^ 
centrations such as cutouts has been of great concern for years in the design 
and analysis of composite structures. Composites with circular cutouts of 
various sizes exhibit the so-called hole-size effects, wtiicli have recently 
received considerable attention from researchers because of their significant 
variations of material strength reduction [1], Thus, the problem of crach exten- 
sion from a hole in composites has been sub:)ected to extensive Investigation 
both analytically and experimentally. Wliile the importance of the problem has 
been well recognized, understanding of the fundamental nature of the problem 
is still very limited due to its Inherent complexities. Various studies have 
been conducted to examine the stress concentration around cutouts in infinite 
and finite-dimensional anisotropic composites and have achieved significant 
success [2-5]. However, research progress in studying crack initiation and 
growth from a hole in anisotropic composite has been relatively slow. The 
difficulties involve not only the complex crack geometry but also the strong 
material anisotropy. In addition, the crack-tip singular behavior and its 
interaction with the hole and the finite-dimensional boundaries of the com- 
posites are inherently associated with the problem. These complications make 
it extremely difficult to be solved mathematically, even for a purely elastic 
case « 

garly work by Bowie [6] has provided solutions for the problem of a crack 
(or cracks) emanating from a circular hole in an isotropic infinite plate sub- 
jected to uniaxial and biaxial stresses. The stress intensity factor coeffi- 
cients, expressible as a function of the ratio of crack length to hole radius, 
are given explicitly in [7] in a tabulated form. Problems of crack extension 


a.** 'iiiM iiB’i iiM II 


2 


from a hole in a flnite-dlmenslonal Isotropic strip have also been invest!** 
gated by many researchers using various numerical and experimental techniques 
[8-10]. However, relatively few studies of the problem for anisotropic com- 
posites have been reported. Waddoups et al [1] were among the first to 
examine the failure of finite-dimensional composites with circular cutouts 
using a fracture mechanics approach. Although Bowie's solution for an iso- 
tropic infinite solid was assumed to be held for the finite-dimensional 
anisotropic composites in the study, results seem to predict the proper trend 
of the fracture behavior of the composites. Recently, a nun^er of alternative 
approaches have been proposed [11-13]. Because of the complex nature of the 
problem and the aforementioned difficulties involved, most of the solutions are 
approximate in nature. ^-Jhile some of the approximate studies provide reasonable 
correlation between solutions and experimental data, no vigorous analysis, to 
the authors' knowledge, of the current problem has been reported. 

This paper describes an analytical method for crack extension from 
a hole in a unidirectional fiber-reinforced composite strip. The unidirectional 
fiber composite, >diich has the simplest laminate configuration, is deliberately 
chosen in this study to provide basic understanding of the problem and to serve 
as a reference for further investigation of more complex multilayered composites. 
General concepts underlying the basic formulation of the current approach are 
discussed in the next section. Solution procedure for numerical computation 
is presented in section 3. The mixed-mode stress Intensity factors character- 
izing the crack-tip behavior are determined explicitly in section 4. Effects 
of fiber orientation on deformation and crack extension are examined for both 
single and double crack cases. Of particular Interest are the energy release rates 
during crack extension; their values are evaluated for several cases to illustrate 
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the fundamental nature of the crack problem. Crack Initiation and growth \ 

from a circular cutout in a multilayered fiber composite give an extremely 
complicated three-dimenaional crack geometry and have to be studied by a 

fracture analysis. Some progress has been made by 
reported elsewhere [ 14 ] * i 
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fully three-dimensional 
the authors and will be 
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2. FORMULATION 


Th© current: approach » using conservation laws of anisotropic elasticity 
and fundamental relationships of anisotropic fracture mechanics » can provide 
rigorous solutions for the mixed-mode fracture parameters without the need of 
solving the corresponding complex boundary value problem* The method has been 
demonstrated to be superior to other approaches in its mathematical rigorous- 
ness and operational simplicity [15], since the conservation integral can be 
evaluated conveniently along a path removed from the crack-tip region. Details 
of the formulation lor the mixed-mode crack problem in anisotropic composites 
have been reported elsewhere [16]; only a brief description of the fundamentals 
of the analysis is given in this section* 

Consider a unidirectional composite plate with rectilinearly aniso- 
tropic elastic properties. Deformation and stress fields in the solid are 
governed by the generalized Hookers law, 


* E a. . 0 . , 

i f* ij j 


i,j - 1,2,6 


wh«e are compllanca coefficients with a^^ - a^^. Compatibility and 
equilibrium conditions of the problem lead to a fourth order governing partial 
differential equation with a characteristic equation of the form 


a^^ - 2a, ^ ^ ( 2 aj ^2 U + a,, * 0 


Lekhnitskli showed [173 that the roots of Eq, 2 are always complex or purely 
Imaginary and occur in conjugate pairs as and p^, p^* for infinitesimal 

deformations of anisotropic elastic bodies, a form of general conservation laws 
for the plane elasticity problem may be written as 
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!h. 


ds » 0 


(3) 


in which S 1b an arbitrary c1ob«4 boundary, enclosing a portion of the continuum 
with unit outward normals n^, and W is the strain energy density. This 
conservation law can be reduced [ig] to the well-known path-independent J- 
integral for an arbitrary path F, which begins on one crack surface, encloses 
the crack tip, and terminates on the opposite face, as shown in Fig, 1, 
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(4a) 


where 




(4b) 


The crack-tip stresses have been proved to possess a classical fracture 

mechanics r"*^ singularity in a linear elastic anisotropic body [19 ) and can 

be characterized by stress Intensity factors K^. and (opening and in-plane 

shearing modes, respectively). For a mixed-mode fracture problem, the correspon- 
ding energy release rates, Gj and G^-j, may be related to Kj and Kjj by 


G,.- 
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Ul ^2 
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and 


G_ 1 * + G-_ 

Total I It 


( 7 ) 
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Using Eqs. 5-7» J can be shown to be teiated to the stress Intensity iactors 
oi a cracked, anisotropic coMposite plate subjected to in-plane loading by 

J » a,, k| + c- . K* K** + K,? (8) 

11 I 1 1 X IX 22 IX ' ' 


“u"- 


*22 ..K + '*2 


tij 


I-(Ui + 


*22 f 1 1 *11 

V-*- V 


It should be noted that Eq. 8 alone does not provide adequate infortnation to 
determine and separately £or the present mixed-mode crack problem. 

If, however, the conservation law for two independent equilibrium states 
[18] is Introduced, additional information may be obtained to remove this 
deficiency. Denoting the two equilibrium states, 1* and 2®, by the super- 
scripts of field variables, the J-lntegral for the superimposed state 0* may 


be expressed by 


where 
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r(k) 




dy - 


d» 


k - 0,1,2 (14) 

(no summation on k) 


u<f» . uf > + u<« 


and is the nutual potential energy deneity, denoted by 

" ‘=ijki “i.j “k.i Saij “i.j \.i 


(15) 


(16) 


Recalling the J-K relationship in Eq, 8 and the orthogonality of the 

ri 2 ) 

stress Intensity factors, M' * can be shown [16 1 to have the form 

K<« K<f + K<« K<i>) + K«) (17) 


The M-lntegral in Eqs. 13 and 17 deals with interaction terms only and can be 
used directly to determine the mixed-mode stress intensity solutions of the 
present problem. 
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3« SOLUTION? PROCEDURE 


Equation 13 together with Eq, 17 provides » in fact* autiicient information 
for determining K<^> end K<« for the current »lx.d-od. enleotroplc frecture 
probiem if proper auxiliary aolutiona are provided. Introducing the first 
auxiliary solution, denoted by the superscript 2a, for the cracked body sub- 
jected to mode I deformation only so that 

j^(2a) . ^ j^(2a) ^ ^ 


Eq.l7 becomes 
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(I, 2a) . ^ ^(1) ^ g(l) 


*11 *'1 


‘12 "n 


(19) 


where the Integral has the form 
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(l»2a) 


r 


. (2a) (1) 

C dv ^ + T^2a) 

'*i,J \,i Ci 9x ^ ^i 


9x 


ds 


( 20 ) 


and In Eq. 20 can be determined easily along a properly chosen 
i 1 

integration path T. Any convenient method such as the commonly used finite- 
element analysis, which can provide accurate far field solutions, may be used, 
u^^®^ and can be calculated from the following relationships: 


(2a) 

“i 


/f ff ^Pi.ui.o) 


( 21 ) 


and 


j(2a) ^(2a) 

*1 ^ij *'j 


( 22 ) 


where 




and f|^^ and g^^^ are given explicitly in Ref. [19], 


( 23 ) 
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Similarly* a seconcl auxiliary aolution* danotcd by the auperacript 2b , for 
pure mode II deformation of the cracked body le Introduced also «o that 


K 


<2b) 


and 


^II ^ 


(24) 


This gives 


» a,„ + 2a«- 

n ^ ‘“22 II 


(25) 


where 


f t 
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(l*2b) 


[ ijkl 1 J k,A ^ 


a,i<2b) a„(l) 

T<1) Zi_ 

^1 3x 1 3x 




(26) 


The quantities and u^^^^ in Eq, 26 may be determined by 

T<2b) - n 

^1 ^IJ "j 


(27) 


whn;re 


and 


.<2b) 


(ID 


- /I 


(28) 


(29) 


( 1 ) 


Thus, Eqs, 19 and 25 provide a system of linear algebraic equations, and K£ 
and can be determined easily. It should be noted that in the process of 
solving for and the integrals, and need to be 

evaluated accurately and explicitly. For a given crack geometry and loading 
condition, this can be achieved by integrating the functionals in Eqs. 20 and 
26 numerically along a convenient path In the far field so as to avoid the 
crack-tip complications. In this paper a finite element analysis is used in 


1 ( 


conjunction with a isecond-order Gaussian quadrature for evaluating the M-int 
in each element (Fig. 2) hy 


N 3 

I l 

n»i n»"l 



and is a weight coefficient in the Gaussian quadrature; 

Cartesian coordinates of Gaussian stations; N, the total number of elements 
associated with the integration path F, and As^ is a segment of V in the n-t 
element . 





XI 


A, RESULTS AND DISCUSSION 

The study of mlxed-niode crack extension from a circular hole in an 
anisotropic fiber composite strip has been conducted successfully using the 
method of analysis mentioned in the previous section. The traction and displace- 
ment fields » T^^^ and involved in Eqs. 20 and 26 for determining 

and may be evaluated along a properly selected path removed from the 

crack tip by a finite-element method* which employs eight-node isoparametric 
elements of the displacement model. Accuracy of current results Is illustrated 
first by examining a crack-hole problem whose solution is available in the 
open literature. Mixed-mode crack-tip stress intensity solutions and the 
associated energy release rates during crack extension are obtained then 
for various cases of interest. 

A. I Accuracy of Solutions 

Accuracy of current solutions for the mixed-mode fracture mechanics 
problem has been studied by examining the problem of a crack emanating from a 
circular hole in an isotropic finite-dimensional plate subjected to a uniaxial 
in-plane stress. Figure 3 shows the geometry and individual dimensions of 
the problem. A crack of length a is assumed to emanate radically at a position (R»<J>) 
on the hole boundary, where R is the radius of the hole and <|> is the direction 
along which the crack extends. For comparative purposes, the problem with 
(|) » 0® was selected for evaluating the accuracy of the solution for the crack- 
hole problem. The analysis, employing only 56 elements in the numerical com- 
putation, gives the result shown in Table 1, which is in excellent agreement 
with the existing solution reported by Stalk and Or ringer [20]. A detailed 
study of the convergence and accuracy of the present method of analysis 


xz 


for mixed-mode fracture problems has been reported elsewliere [15,16!] . The results 
indicate that, in general, solution accuracy within one percent deviation from 
the reference ones obtainable in literature can be achieved easily using only 
one-third of elements in other highly sophisticated, singular finite-element 
analyses [20,21]. The accuracy of the current solutions is further reported 
to be insensitive to the degree of discretization in the numerical calculation. 
This situation is expected, since the present approach bypasses the highly 
sensitive and complicated crack-tip region and evaluates the conservation 
Integral along a selected path in the far field, where field variables are 
generally smooth. 

4.2 Cracks Emanating from a Circular Hole in Composites with Various Fiber 

Orientations 

More complex fracture problems of a crack (or cracks) emanating from a 
circular hole in off-axis unidirectional fiber composite strips are examined 
in this section. Since the response of the problem is antisymmetric in nature, 
two cracks of equal length, a, are assumed to be initiated from the boundary of 
the hole in the composite, where high stress concentration occurs. The corre- 
sponding single crack case is also studied for comparative purposes. Figure 3 
shows the geometry of the cracks and the circular hole located at the center of 
a finite-dimensional 2b x 2L anisotropic panel subjected to uniaxial nominal 
stress, 0 ^, along the y-direction. The material axes, denoted by 1 and 2, of 
the plate are oriented with respect to the structural axes x and y by an angle 6, 
as indicated in the figure. The local stress field around the hole (without a 
crack) in the finite-dimensional cr^osite strip is solved first to determine 
accurately local stress concentrations by an anisotropic finite -element analysis 
developed in the present study. The Tsai-Htll stress failure criterion [22] is 
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Introduced then into the analysis to determine the locations, (R,^) and (R,7f+^), 
of crack initiation from the edge of the hole by 

■ [in . ■ ~w iTj + -s2" “ ^ 

where subscripts 1 and 2 refer to the stresses in the material coordinates, 
and X, Y, and S are the longitudinal, transverse, and shear strength of the 

unidirectional composite, respectively. While other criteria can be Imple- | 

I 

mented easily in the current analysis, no attempt has been made to evaluate | 

the accuracy of different failure theories. For an E-glass/epoxy composite | 

strip panel, X, Y, and S have the following values [23]: | 

X * 150 X 10^ psl 

‘i 

I 

Y ® 4 X 103 psl ? 

• i 

I 

!■ 

and S * 6 X 10^ psi | 

Values of (p for various degrees of material anisotropy (l.e. fiber orientation) I 

I 

of the composite are obtained by Eq. 31 in the present analysis. Figure 4 | 

gives the location of crack emanation as a function of 6 for the E-glass/ 
epoxy composite indicating that <P decreases appreciably with Increasing 6. 

I 

It is a common practice in studying fracture of the unidirectional composite | 

to assume that cracks extend, after Initiation, along a direction parallel to | 

the major material axis, l.e. the fiber direction, since the trans-fiber 

strength of the composite is generally much lower than that along the fiber | 

I 

direction. Now consider the problem of two cracks of equal length, a/R *= 0.4, j! 

I 

emanating at (R,(J>) and (R,u+(|>) from the edge of the hole in the glass/epoxy | 

I 

composite, where <|i depends upon the fiber orientation as shown in Fig. 4. | 

Fracture parameters such as mixed-mode stress intensity factors, Kj. and | 

I 

I 

I 

\ 
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which describe the crsck-tip deformation and fracture of the composite » are 
calculated explicitly for cases of various fiber orientations t For compara- 
tive purposes, corresponding single-crack problems are also analyzed. Results, 
expressed in terms of Kj. / (o^/iTa’) and both single and double 

crack cases, are reported in Fig* 5. The associated energy release rates G 
are determined also in Fig. 6. The results indicate that K| Increases contin- 
uously with the increase of fiber orientation and that rises to a maximum 
as 6 approaches 45^ and then gradually decreases to zero. The opening mode 
stress intensity factors determined for most cases are generally larger than 
that of the shearing mode. For the particular crack length studied, differ- 


ences in K^, Kj.^, and 6 between the single and double crack cases do not seem 


to be appreciable. Nonetheless, a significant discrepancy may be expected if 


the a/R ratio becomes larger, which will be discussed in the next section. One 
extreme situation is the case of 0 90^ in which the crack cs tension is shown 


to occur along the direction normal to the loading, i.e. (J) « 0^, as expected. 


This gives a purely opening mode fracture,!, e. 0.00, due to the symmetry 

of loading and geometric conditions* 


4*4 Crack Extension in Off-Axis Unidirectional Composites 

Problems concerning crack extension from the edge of a hole in the com- 
posites are also examined in detail in this study. The case of an off-axis 
unidirectional glass/epoxy with fibers oriented along 0 » 45° is presented 
here for an illustrative purpose. The geometry and loading conditions are taken 
the same as previous examples (Fig. 3) . Both single and double crack cases are 
considered again* Solutions for the problem have been obtained by an incre- 
mental scheme using the method of analysis described in section 2. Results 
given in Figs. 7 and 8 illustrate the fundamental behavior of crack extension 
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from a hole in the off-axis composites. Effects of the crack length on stress 
intensity solutions and energy release rates are shown clearly in the figures. 
Values of Kj./(v^ aj) remain relatively constant for a wide range of the a/R 
ratio studied and change significantly for both smaller and larger cracks in 
the single crack case. For the double crack case, K^/(/if 0 ^) increases mono- 
tonically during crack extension. It is interesting to note that, for both 
single and double crack cases, varies approximately linearly with 

the crack length. Slight deviations from linearity in 0 ^) are observed 

as the crack (or cracks) extends beyond a/R » 1. In general, the results indi- 
cate that all Kj., and G have larger values in a double crack case than tliose 

in the corresponding single crack problem, presumably due to the larger ’’effective'* 
crack length. Moreover, the larger the crack length, the greater the difference 
between the two cases. Composites with other fiber oriencations may 

exhibit different values of K’s and G*s during crack extension, but fundamental 
features of the crack-tip deformation and fracture behavior are essentially 
the same. 

Extension of this method to more complicated fracture problems of composites, 
such as viscoelastic anisotropic crack analysis, can be easily achieved. These 
further developments have been very successful and are reported elsewhere [24]. 
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5» SUMMARY AND CONCLUSION 

An analytical investigation of cracks emanating from a circular hole in 
an off-axis unidirectional fiber-reinforced composite has been conducted. The 
method of analysis employed in the study is formulated by using conservation 
laws of elasticity and fundamental relationships in anisotropic fracture 
mechanics. The current approach provides a very convenient and accurate means 
to examine the complicated crack behavior* when used in conjunction with a 
suitable numerical scheme such as the finite element method. The problem is 
eventually reduced to a system of linear algebraic equations in mixed-mode 
stress intensity factors. Fracture parameters* describing crack-tip defor- 
mation and fracture in the composite* are obtained explicitly. Effects of 
material anisotropy and ct'ack/hole geometry are examined also. Of particular 
interest are the energy release rates associated with crack extension; their 
values are evlauated for various cases. 

The results obtained in this study lead to the following conclusions: 

(1) The current method of analysis is demonstrated to be superior to other 
approaches in its operational simplicity and solution accuracy for In- 
vestigating mixed-mode crack problems with extremely complex crack hole 
geometry and material anisotropy. 

(2) Mixed-mode stress Intensity factors and energy release rates associated 
with the cracks emanating from a hole change very appreciably with fiber 
orientation in the composite. K^, and G increase monotonically with 
increasing 0; but reaches its maximum at 0 ■ A5°, and then decreases 
gradually as 0 increases further. 
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(3) During crack extension from a circular hole in the composite with a 
given fiber orientation, the opening mode stress intensity factors are 
found to be relatively constant for a certain range of the crack length (i.e,, 
a/R < 1»0) in the single crack case, but increase rather rapidly In the 
double crack case* increases linearly with the crack length, but 
deviates gradually from linearity as the length of the crack exceeds the 


hole radius, l.e a/R > 1, for both cases. 

(4) The energy release rates and mixed-mode stress intensity solutions for 
single and double crack cases do not have much difference for problems 
with small crack-length- to^hole-radius ratio i.e. a/R < 0.5* However, 
the difference increases and becomes very appreciable as a/R increases 


during crack extension, 


I 
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TABLE I 

STRESS INTEHSITV FACTORS FOR A CRACK EMANATING FROM A CIRCULAR HOLE 
IN AN ISOTROPIC FINITE-DIMENSIONAL PANEL* 




*^XI 

Stalk and Orringer [20] 

i.m 

0.000 

Present Solution 

2,333 

0,000 


*a/R » C.4* 0®, 56 eight-node isoparametric oXements 


1 


n 


LISf OF FIG0RE CAPTIONS 


Figure I» Goordinatee and Path for J^^IntegraX 

Figure 2, Gauaaian $tatlone, Finite Elenent Heah» and Path for M-Integral 
around Crack Tip 

Figure 3« Cracka Enanating from a Circular Hole in an Off-Axis Unidirectional 
Fiber Composite with Material Axes Oriented along 1-2 Directions 

Figure 4* Location of Crack Emanation from a Hole in Glase/Epoxy Composites 
with Various Fiber Orientations, a/R«0.4 


Figure 5«> Mixed-Mode Stress Intensity Factors of Cracks Emanating from a 
Circular Hole in Unidirectional Glass/Epoxy Composites with 
Various Fiber Orientations, a/R«0.4 

Figure 6, Energy Release Rate of Cracks Emanating from a Circular Hole in 
Unidirectional Glass/Epoxy Composites of Various Fiber Orien- 
tations, a/R-0.4 


Figure 7* Mixed-Mode Stress Intensity Factors vs. Crack Length in Unidi- 
rectional Glass/Epoxy Composites with 0 •• 45 " 

Figure 8. Energy Release Rate vs. Crack Length in Unidirectional Glass/ 
Epoxy Composites with 8 45® 




FIG. 1 


COORDINATES OF THE CRACK AND PATH FOR 0 INTEGRAL. 
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16. 2 GAUSSIAN STATIONS, FINITE ELEMENT MESH AND 
PATH FOR M INTEGRAL AROUND CRACK TIP. 
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LOCATION OF CRACK EMANATION FROM A CIRCULAR HOLE IN 6 
COMPOSITES WITH VARIOUS FIBER ORIENTATIONS. a/R « 0.4 
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FIG. 5 MIXED-MODE STRESS-INTENSITY FACTORS OF CRACK(S) 
EMANATING FROM A CIRCULAR HOLE IN UNIDIRECTIONAL 
GLASS/EPOXY COMPOSITES WITH VARIOUS FIBER ORIEN- 
TATIONS, a/R = 0.4. 
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FIG. 6 ENERGY RELEASE RATES OF CRACK(S) EMANATING FROM A 

CIRCULAR HOLE IN UNIDIRECTIONAL GIASS/EPOXY COMPOSITES 
WITH VARIOUS FIBER ORIENTATION, a/R = 0.4. 
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FIG. 7 MIXED-MODE STRESS INTENSITY FACTORS VS. CRACK 
LENGTH IN UNIDIRECTIONAL GLASS/EPOXY COMPOSITE 
WITH 9 = 45*. 
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FIG. 8 ENERGY RELEASE RATES VS. CRACK LENGTH IN 
UNIDIRECTIONAL GLASS/EPOXY COMPOSITE WITH 
0 - 45\ 



